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ABSTRACT 


Accurate charting of marine fog distribution is made 
feerreult because of the sparseness of transient-ship obser- 
mavens. Ihe meteorological satellite offers the possibility 
of improving fog diagnosis by specifying areas (vice a dis- 
crete number of points) of marine fog. 

mats study explores the feasibility of using meteoro- 
logical satellite data, visual and/or infrared, as means of 
discerning the presence of marine fog. Transient-ship and 
satellite data from a summer-season period in 1973, eastern 
North Pacific Ocean, served as data base. A qualitative 
comparison is made between NOAA-2 satellite imagery and 
verifying marine fog, as determined from the synoptic-time 
Ship data; qualitative guidelines for diagnosing the presence 
of marine fog from the visual and infrared imagery were 
Hornivbated. A quantitative determination of the presence of 
Pomme LOC. Using yYisual and infrared NOAA-2 digital satel- 
lite data, is presented, from which a technique is developed 


Puno vectively Giagnosing regions of marine fog. 
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I. INTRODUCTION 


A. BACKGROUND 

MAmGincerOO BUS an ever-present menace to both commercial 
and military shipping. The ship master and ship owner are 
threatened with costly schedule delays and losses of property 
and personnel because of the possibility of collision during 
periods of low visibiiity in marine fog. Likewise, military 
mival £orces afloat and naval air forces are hampered by fog 
iomamc evolutions requiring unrestricted or near-unrestricted 
Visibility. For example, aircraft carrier flight operations 
(launch and recovery), antisubmarine warfare maneuvers (ship 
Based and shore based), tactical air-strike operations, and 
underway replenishment, all require unimpaired horizontal 
Visibility. Wheeler (1974) has documented the importance of 
marine fog to naval operations historically and has submitted 
representative statistics of losses in dollars and lives 
suffered by the United States Navy as a direct consequence 
See low visibility due to marine fog. 

Realizing the importance of marine fog to maritime ship- 
ping, various Department of Defense (DOD) activities have 
cieaccdunii rescarch imnvolyine the analysis, modeling, and 
forecasting of fog over open oceans. Englebretson (1974) 
ists the participating DOD activities and outlines their 
respective research goals. The Naval Postgraduate School 


Se oyecroup, inesparticular, is primarily focusing on the 
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analysis and forecasting of marine fog. A portion of the 
group has addressed the problem of marine fog climatology 
emoedevetoped a method of synthesizing ships't surface synop- 
Mmemreports for the purpose of assessing the climatological 
frequencies of marine fog (Renard, Englebretson, and 
Daughenbaugh, 1975). The scheme continues under development. 
eiimatologyers generally dependent on large volumes of 
data, hopefully near-uniformly distributed both spatially 
and temporally. However, a climatology of any given para- 
meter over regions where the data coverage ranges from sparse 
to relatively abundant has certain obvious biases that may 
render misleading results. For example, for a geographical 
Pi@eawoucside of the normal shipping lanes where few or no 
ships pass, a general fog climatology might show the area 
"fog-free."' In reaiity, the climatologist is probably not 
eens there are no occurrences of fog in that area, but 
merely that there are no "reported" occurrences of fog in that 
Parone Otnicr proplems innate to climatologies compiled 
PeommeamwenctTONnal Synoptic ship data are described by Renard, 
Pmerebretson, and Daughenbaugh (1975). However, it should be 
pointed out that climatologies thus derived are probably the 
meoewacnievable for the current “state of the art." 
PoateeiemmetcOrolocist really desires then 1S some 
sort of data-taking apparatus that would give better temporal 
and spatial continuity to an observable phenomenon. Such an 
instrument might be the meteorological satellite provided the 


phenomenon being observed is, in fact, remotely discernible 


ual 





Bpesaecellite sensors. Several satellite meteorologists in 
Bitcmracid have suggested the satellite approach for detecting 
marine fog while, for example, Anderson, et al (1969), who 
Mreat the subject of weather-satellite applications in depth, 
memnot even address the subject of open-ocean marine fog. 

If it can be determined that marine fog is discernible by 
Zeereeec SeUuSOnS, Complete real-time description of fog 
Mest ributcion may at last be a reality. Hence, the time has 
mone In the marine fog studies to determine whether or not 
marine fog is observable by current sensors on board 


meteorological satellites. 


B. OBJECTIVE AND APPROACHES 

The objective of this study was to determine the feasi- 
bility of using meteorological satellite data, visual and/or 
infrared in any format, as means of discerning marine fog, 
Pmemaimcensity, and its spatial extent. There were basically 
two approaches taken to achieve this objective. The first 
was a qualitative and semi-quantitative interpretation of 
Saucllytec imagery in order to derive and set forth qualitea- 
tive guidelines by which trained meteorologists at field 
activities could establish the existence and extent of marine 
fog from visual and/or infrared satellite imagery with some 
degree of positive skill. The second was the processing and 
Wiieteakive interpretation of digital satellite data in 


order to establish the importance of applying numerical 


EZ 





techniques to satellite-derived data as a means of more 
accurately and objectively outlining regions where marine 


fog exists. 


C. DATA 
imeeoynopcic Ship Data 

Mice ne shitpesyNopeleatime weather ODSErVations, 
macembastcally as ‘'sround truth'’ for the existence or non- 
existence of marine fog, were obtained from two sources. 
The first source was the National Climatic Center (NCC), 
Asheville, North Carolina; the data were provided by the 
Naval Weather Service Detachment (NWSD), Asheville, North 
Carolina, on magnetic tape in a format known as Tape Data 
Mamiiy-11 (TDF-11), Surface Marine Observations. The second 
source was the Fleet Numerical Weather Central (FNWC), 
Memterey, California; these archived observations were also 
Waevided On Magnetic tape in a format referred to as SPOT. 

it Piewas provided for summer (June through Septem- 
Cero s, ollce the data were So nearly current at the 
mianenOr Lrequest, the usual NCC quality checks had not been 
made. This absence of quality control probably rendered the 
Gitasmonke Trepresentative OL “Uunpurged,' real-time synoptic 
Sesctrvations. Computer sottware was written to extract 
pertinent information from the 140 character TDF-11 format 
oma rerormat it in standard Ship synoptic weather code. 

SPOT was pre-processed at FNWC to provide a printed 


copy of synoptic ship observations in an abbreviated, easily 
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readable format. These synoptic data were randomly checked 
against the TDF-11 to ensure that one source was not collect- 
ing more or different data than the other. It was found 
that, generally, for (militarily) unclassified data the two 
mets were nearly identical. The eventual purpose of this 
Somparison was to determine the aptness of synoptic charts 
meotted from the FNWC SPOT using the FNWC CDC 6500 digital 
meeurer and on-line Varian electrostatic plotter. 
eo ce eves Pressure Analyses 
National Meteorological Center (NMC) Northern 
Hemispheric sea-level pressure analyses on microfilm were 
Peoavaded by the Environmental Prediction Research Facility 
(EPRF}, Monterey, California. The NMC analyses were used 
for determining general synoptic weather patterns, particu- 
tardy major cyclones, frontal activity, and surface wind 
flow. In addition, the plotted synoptic ship observations 
were employed as a tertiary source of raw synoptic data. 
feeeececorological Satellite Data 
Pe NOl eee ulOs- bis Satellite Data 

NOAA-2 data were obtained from two sources. EPRF 
Puovidcasc. <0 photographic enlargements of hemispheric 
Gieeitally composited mosaics for the scanning radiometer 
daytime visual (SRVIS), daytime infrared (SRIR-DAY), and 
nighttime infrared (SRIR-NGT) observations. Composited 
digital SRVIS and SRIR data on magnetic tape were provided 
by the NCC through the NWSD, Asheville, North Carolina. Both 


the mosaics and the digital data were originally processed by 
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mere National Environmental Satellite Service (NESS), Washing- 

mone GCG. NOAAN=2Z data were selected for this study primarily 
because of the worldwide availability of NOAA satellite read- 

meee (direct or archived) information to civilian and military 

activities. 

The NOAA-2 satellite is in a near-circular early 
morning/early evening orbit at an average distance of 783 nmi 
(1451 km) from the earth. Equatorial satellite sub-point 
crossing times occur near 2100 local northbound and near 0900 
hoeate southbound, The orbital period is 115 minutes or 
roughly two hours. 

The NOAA-2 scanning radiometer has a visual 
channel with a spectral response of 0.5um to 0.7um and an 
infrared window channel with a spectral response of 10.5um 
mote eoum. Lhe SRVIS channel has a 2 nmi resolution at the 
nadir point and the bandwidth of the spectral response is 
iMerre ly cCOMmtinuous Within the visual spectra. The SRIR 
Channel has a 4 nmi resolution at the nadir point and is a 
relatively “clean" atmospheric window, oriented to the long 
wavelength side of the 9.6um ozone absorption band and to 
the short wavelength side of the 15um carbon dioxide 
absorption band (Cogan and Willand, 1974). Some absorption 
by water vapor does occur toward the long wavelength end of 
the IR spectral interval. (See Figure 1) The usefulness of 
ngs water vapor Sensitivity, especially in NOAA series 


Very High Resolution Radiometer (VHRR) IR data, for the 





detection of lower tropospheric moisture and areas of 
possible fog formation, has been postulated (National Weather 
emyvuce, Western Region Headquarters, 1974). It should be 
noted, however, that, although the absorption due to water 
Vapor in the higher wavelengths of the NOAA IR band is sub- 
Stantial, it is still considerably less than the absorption 
for the atmosphere as a whole, (See Figure 1) 

iiemacdelky sNOAA=Z mesartes used in) this study were 
Sebsaiialby prepared by high speed digital computers. SRVIS 
pdeoktR data Obtained during the daylight portion of the 
Srbiemti.e., southbound) are used in the O0900L mosaics; SRIR 
tetas oObtained during the nighttime portion {1.e., northbound) 
are used in the 2100L mosaics. The scan-line analog signals 
from the satellite are digitized using analog-to-digital 
conversion equipment. The digital information is then 
processed, earth located, and overlaid on a standard polar 
Stereosraphic projection. Where data overlaps on consecutive 
Orbits, the more current data are used. <A more complete 
description of the content and preparation of the mosaics 1s 
Sontained In any volume of Key to Meteorological Records 
Documentation (KMRD) No. 5.4 (NESS, 19XX). 

The digital NOAA-2 data were provided on magnetic 
Mapes. Lach tape contained either the SRVIS and SRIR-DAY 
or SRIR-NGT for one date. The SRVIS/SRIR-DAY data were 
configured in a 2048 x 2048 mapped array for the Northern 
Hemisphere; documentation for each tape was contained in the 


Pict VecemanOr the data File; and records Z through 2049 
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Someained one line each of the 2048 count values, ranging 
meemeoetoO 255 (0 to 377 octal). The SRIR-NGT data were 
provided for the Northern and Southern Hemispheres; the 
first file contained the documentation record; the second 
file contained 2048 records of the 2048 x 2048 Northern 
Hemisphere SRIR-NGT count data; and the third file contained 
2048 records of Southern Hemisphere data. A schematic of 
the tape layout for both sets of tapes is provided in Figure 
2. For SRVIS, low order values indicate dark and high order 
values (254 maximum) indicate bright; for SRIR, low order 
values indicate cold and high order values (254 maximum) 
mreteate warm. Value 255 (377) is reserved for denoting 
"void" areas which appear either as the peripheral back- 
meoundeor 2S data-void regions. At NESS, current digital 
Weocessing Capabilities include such things as: corrections 
for atmospheric absorption (limb-darkening), conversion of 
IR radiances to equivalent blackbody temperatures, rectifi- 
cation to different map scales, and gray shade enhancement 
for IR imagery. <A more detailed discussion of these and 
Giemeneprocessing techniques at NESS are discussed by Anderson, 
et al (1969). 

EPRF has developed equipment and support computer 
software for a similar type processing of NOAA series data at 
Naval Weather Service Fleet Weather Centrals which possess 
Cee oAVONpeocessing capabilities. The Direct Readout Data 
Comvertenr (DRDC)) as it is referred to, and the DRDC support 


software package are designed to digitize, process, and 
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@oplay the NOAA series DRSR data. One particular processing 
memrime, Called “Slice,” produces an output picture which 
contains only a few discrete gray shades, each related to a 
mecitic temperature range (IR only) or brightness range 
(visual only). It has been speculated that this processing 
mouwerme COUuld be useful in marine fog studies. A complete 
documentation of the DRDC and DRDC software is provided by 
momolt (1974). 
Dt oumromnenr ly NAPPiwoatellite Data 

Sate wmtewadtaerrone che sNerense Meteorological 
eaeellite Program (DMSP), formerly referred to as Data 
Acquisition and Processing Program (DAPP), were chosen for 
Study because of the planned widespread availability of 
DMSP satellite direct-readout at Department of Defense 
weather activities both ashore and afloat. In addition, DMSP 
communications and ground processing systems are configured 
so that a usable product is available within five minutes of 
the termination of the transmitted data stream (Meyer, 1973). 

DMSP satellites operate from an altitude of 450 
nautical miles in sun synchronous orbits. Two satellites 
are normally maintained in orbit, one in the noon-midnight 
Paceemandethe Other in the carly morning-early evening orbit; 
the DMSP sensors are designed to function at any local sun 
Mice ine satellites cross the equator northbound at 
approximately 0830L/1230L and southbound at 2030L/0030L. 
The DMSP satellite visual sensors have a spectral response 


in the region between 0.4um to 1.lum which covers the 
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near-ultraviolet to near-infrared portion of the solar 
Spectrum; this region was selected to optimize the distinc- 
tion among clouds, ground, and water (Meyer, 1973). One of 
the two visual sensors has a 2 nmi resolution (H-format) for 
global coverage and the other sensor has 1/3 nmi resolution 
(V-format) for limited area coverage. In the preliminary 
femaies for this project, the unique DMSP visual imagery 
Showed promise as a tool for detecting lower tropospheric 
moisture, and/or aerosols and condensation nuclei; this 
eet 1S discussed in more detail under "Procedures." 

The infrared sensors cover the broad bandwidth 
of 8.Qum to 13.Q0um which spans the so-called "atmospheric 
window."' However, it is essentially a "dirty" window, being 
"contaminated" by water vapor, ozone, carbon dioxide and 
femesols (Cogan, Tripp, and Isaacs, 1974). (See Figure 1.) 
The usefulness of such a broadband IR spectral response, 
especially as "dirty windows" for lower tropospheric moisture 
detection, has been proposed (Cogan and Willand, 1974; NWS, 
Western Region Headquarters, 1974). Like the visual sensors, 
one IR sensor has a 2 nmi resolution (I-format) for global 
coverage and another IR sensor has 1/3 nmi resolution (W-format) 
for regional coverage. 

Unlike the NOAA series satellites which transmit 
Seisca iitLormation as an analog signal, the DMSP satellite 
data are digitized onboard and transmitted in digital format 
PidenMmicgwane data digitally processable upon receipt. For 


example, the IR output voltages are converted onboard the 


L3 





Zieeiite toad linear function of the equivalent blackbody 
Memperature of the scene (to give a uniform temperature 
resolution) within a temperature range of 210-310K and eens 
transmitted to earth for direct processing (Blankenship and 
eavace, 1974). In addition to the ability to process IR 
data at various temperature ranges and thermal resolutions, 
the current DMSP processing and display equipment has the 
capability of so-called "threshold processing" whereby switch 
selectable threshold values may be established at three IR 
memperatures within the 210-310K range. The output image 
contains four shades of gray corresponding to temperature 
intervals within the chosen interval. (This is very similar 
eemePRe's "Stice'' technique described on page 18.) It has 
been speculated that this threshold processing might be 
extremely helpful in lower tropospheric moisture and fog 
Studies. Blankenship and Savage (1974) and Meyer (1973) 
Gece a rather brief but concise overview of the DMSP sensors 
and electro-optical processing capabilities. 

Dismedata ton stne North American and the North 
Ltetrre Ocean areas for June through September 1973 were 
obtained from the Air Force Global Weather Central, Offutt 
Air Force Base, Omaha, Nebraska, petore they were system- 
ically degaussed. Also obtained were the accompanying 
ungridded image transparencies of the H- and I-configured 
data. The AFGWC converted the raw encrypted H and I DMSP 
Tuciluce datasto mapped data format (enhanced and rectified) 


with accompanying photographic imagery of the gridded 


20 





Spin-scan display on a 1:15 million polar stereographic pro- 
jection. The conversion of the raw encrypted data to mapped 
data format renders the data computer processable and earth 
Mocatable. 


c. NIMBUS-5 Temperature-Humidity Infrared Radiometer 
(THIR) Data 


NIMBUS-5 THIR gridded film strip imagery for the 
Oem (6.5-7.0um) and 11.5um (10.5-12.5um) channels were 
obtained for EPRF. The 6.7um water vapor channel gives infor- 
mation on the moisture content of the upper troposphere and 
stratosphere and thereby jet stream and frontal activity 
tiecdanrdmopace Flight Center, 19/72). For this reason, the 
6.7um channel imagery was only briefly perused and then dis- 
regarded for this study. The THIR 11.5um window channel 
povided day and night cloud top or surface temperatures (for 
clear column radiances). The 11.5um channel has essentially 
the same spectral response as the NOAA-2 SRIR sensor, 1i1.e., 
miesoum to 12.5um. According to Cogan et al (1974), the 
NOAA-2 SR bane or Pim Demtch menanecidt or the NIMBUS-5 
fiekeliesum channel im terms of less atmospheric attenuation. 
Even though atmospheric attenuation, especially water vapor 
Spsoupeion, Might be beneficial in fog detection studies, the 
difficulty in viewing the gridded film strip imagery precluded 
its use in this study for other than as a verifying aid in 
questionable areas. Moreover, the NIMBUS series satellites are 
not designed for direct readout capabilities and thus cannot be 


Escderone real-time meteorological operations. 
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Figure 1 has been provided for a comparison among 
the spectral bands covered by the NOAA-2, DMSP satellites, 
pom NIMBUS-5 THER Il.5um sensors relative to the atmospheric 
meeOrption spectra. Figures 3 and 4 are samples of the 


THIR 11.5um channel, NOAA-2 SRVIS and SRIR-DAY, and DMSP 


weet, and I imageries. 
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eee ROCEUURES 


fee oOLELECTING THE GEOGRAPHICAL REGION AND TIME PERIOD 

Pies geoeuraphnical area for this study was chosen to be the | 
North Pacific Ocean from 30N to 60N latitude between 115W 
and 180 longitude. (See Figure 5.) This region was selected 
memncarly GCO1ncide with part of the area that Renard, 
Englebretson and Daughenbaugh (1975) used in develcping 
climatological marine fog frequencies. This North Pacific 
Ocean area afforded a suitable domain in which both maxima 
and minima of fog frequencies occur, as shown in Figure 6. 
Choosing the area of study in this manner insured that the 
Miield-of-view'' would be large enough to include the critical 
boundaries between fog and no-fog. 

The time period of this study was selected based on a 
pre-screening of National Meteorological Center (NMC) six- 
hourly synoptic surface charts (i.e., 0000 GMT, 0600 GMT, 
ied Gits  hs00 GMT) for the months of June through September 
Povo. Periodsewhen £62 existed over a significant portion 
of the study area, as revealed by the transient-ship reports 
and/or synoptic features, were then noted. These selected 
periods were compared to the dates and times for which data 
existed simultaneously from DMSP, NOAA-2, and NIMBUS-5S 
Tie eet tCsmliiinerepard te the longest continuous period 


Baomrirechmedatasexisted trom all Satellite sources, the Limited 
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period of 17 July through 6 August 1973 was selected for 
study. 

iitest: ls iitinon polar stereographic projection was 
chosen as most practical for displaying the NMC sea-level 


pressure analyses, the NOAA-2 mosaics, and the DMSP imagery. 


Bee ERCCESSING TRANSIENT-SHIP SYNOPTIC-TIME DATA 
me OnLeing and Categorizing IDF-1] Data 

Once the [PF-ll data had been reconstructed into the 
Poenvenient standard synoptic format for ships, only reports 
within the time-frame of the study were retained and printed. 
Additionally, for easier use, the ship reports were sorted 
imo niinary SynMoptic reporting times, 1.¢e€., 0000 GMT, 0600 GMT, 
1200 GMT, and 1800 GMT. Data for other time periods were 
preserved separately in order to provide temporal continuity. 
The total number of reports for the primary reporting times 
Metaied 5053. 

Computer software was then written to scan the ship 
reports within the primary synoptic periods for fog or 
fog-related parameters, categorize them according to Table I, 
idee tiem am the format displayed in Figure 7. This 
codification was performed to provide overall fog statistics 
[OLE iecm send yeperrod sand 2350 to aid in selecting a shorter 
time period within the 21-day primary period for a more in- 
tepenectldy ss lhe Sortie of the foe reports revealed that 
present-weather fog reports comprised 20.25% of the 5053 


reporting ships. More generally, the categorization of the 
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fog reports revealed some remarkable internal inconsistencies 
within the synoptic reports which reflected upon the capa- 
bility of establishing the existence of cloud layers above 
meported fog and hence any higher cloud "contamination" that _ 
might be expected in satellite data. For example, a common 
problem occurs in the reporting of total cloud amount (N) 
consistent with the fog-type which is reported. In numerous 
Gases, code figure 8 (i.e., 10/10 cloud cover), vice 9 
(i.e., sky obscured), is reported for "N'" while code figure 
tomeeo, 47, Or 49 (fog, sky not discernible) for "ww" is 
Smiitaneously reported. Also, there are cases where code 
meoure, 9 1S reported for "N'" and clouds at one or more 
levels are reported (i.e., for Cy; Ci: C,)-. In any event, 
Pomona Ships reporting fog also reported an cbscuration 
of the sky or a type of fog that implied that the sky was 
obscured. The task of establishing some "ground truth" of 
higher level clouds (types and amounts) with the existence 
promt sea (i.e... Stratus within 33 feet of the surface 
eee water) 1S Virtually impossible. Thus, the hope of 
establishing cloud types a priori from transient-ship 
weather reports for satellite observations was abandoned as 
MeromenicnDnhGccdure dmethis study. However, this did not 
prove to be restrictive except when fog existed and was 
reported within frontal bands where multilevel-clouds exist. 
2. Plotting Synoptic Data from SPOT 
CanpulitcumcOrtwane was written for the FNWC CDC 6500 


POmolount nosis weather data for the 84 six-hourly, 
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People uciartsscovering the 2l-day study period. The soft- 
ware was designed to extract SPOT data from a FNWC climato- 
logy tape and display the data on the appropriately 
meomeraceae!: 15 million polar stereographic background. The 
data were color-coded to depict present-weather fog, past- 
weather fog, or cloud-types (for ships not reporting fog) 
and specific areas of fog and no-fog were subjectively 


generalized from the synoptic reports. 


fee REPARING SEA=LEVEL PRESSURE ANALYSES 

iemnite sea-level pressure analyses were provided on 
microfilm. In order to more readily use the charts, the 
microfilmed analyses were photographically enlarged to 1:15 
MeeirOne s bY thus enlarging the analyses, plotted synoptic 
data and/or satellite imagery could be easily compared to 
various features and parameters in the analyses. The fog/no- 
mOoOwareas were SUbjectively adjusted, as necessary, to 
wmopetcurecatUures Particularly near fronts, along the per- 
iphery of anticyclones, and in areas under southerly flow 
across the cold waters near the Aleutian (subarctic) ocean 
Gumien t. 
eee ROCESoING NOAA-2 MOSAICS AND THE Bebe Or wNOAA- 7 

IMAGERY FOR DISCERNING MARINE FOG 

In order to view the NOAA-2 mosaics on the same scale as 
the plotted transient-ship weather data and sea-level pressure 
analyses, it was decided to photographically enlarge each of 


thewokVilo., GRIR-DAY> and SRIR-NGI mosaics to 1:15 million 
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ieee suuldy area. There are approximately 99 digital 

Sean lines per centimeter advancing from the top to the bot- 
tom of the normal 20.7 x 20.7 cm mosaic. Since the enlarge- 
ment process required a 7.8 expansion of the mosaic, the 
enlarged imagery has roughly 13 scan lines per centimeter. 
As a result, the photographically enlarged imagery showed 
some degradation in the sharpness of the gray shades. 

There were a few cases where this degradation was aggrevated 
apparently by other processing problems including the initial 
mosaic reproduction. The result was that some of the 
reproductions were generally too dark within the study area. 
This factor was influential in the subsequent choice of the 
Narrower time period for a detailed study. 

In the area of concern, i.e., 30N to 60N between 115W 
ama iG0, there are primarily four NOAA-2 passes which are 
eventually digitized and combined to create the day compos- 
ites (SRVIS and SRIR-DAY, satellite southbound) with 
approximate ascending nodal times of 1648 GMT, 1848 GMT, - 
2048 GMI, and 2248 GMI (stepping east to west on successive 
orbits). This means that the time when the satellite has a 
subpoint of 60N is roughly 40 minutes after ascending node 
(AAN) or 1728 GMT, 1928 GMT, 2128 GMT, and 2328 GMT. The 
Saucelrtre tien Crosses SUN approximately ten minutes later. 
mpulee the westernmost pass (AN 2248 GMT) is used only in 
pilemextrche northwest portion Of the area of study, the 
primary synoptic data used for the day mosaics were for 


1800 GMT with the 0000 GMT data being used secondarily for 
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Sompimurcy. | oimilarly for the night composites (SRIR-NGT, 
satellite northbound), the times of ascending node were 
0448 GMT, 0648 GMT, 0848 GMT, and 1048 GMT. The satellite 
subpoint crosses 30N ten minutes later and crosses 60N at 
approximately 0508 GMT, 0708 GMT, 0908 GMT and 1108 GMT 
(stepping east to west on successive orbits). Thus iene 
0600 GMT and 1200 GMT synoptic data were the primary and 
Secondary verifying "ground truth", respectively, for the 
night composites. 

A comparison was then made between the verifying 
synoptic data and the satellite imagery. Gray shade vari- 
melons eit the Satellite imagery were noted in areas where 
fog was reported particularly along the discernible boundaries 
of varying gray shades which coincided with fog/no-fog bound- 
Penccmiiterrea rom the verifying synoptic data. Each of 
Eiemiewenty-one Simultaneous SRVIS and SRIR-DAY and the twenty- 
one SRIR-NGT enlarged mosaics were compared to the verifying 
ship data and significant features were documented. At the 
Same time, a subjective assessment was made as to the overall 
iucderry oF each enlarged mosaic and to the quantity of 
verifying data; this information was used as an additional 
Times hcCectane a5Sub-Study period from the 21-day period. 

Generally, it was found that, for areas unobscured by 
iuaihemancGmnrene leveleclouds, at was possible tomick out 
areas of low level clouds and further distinguish stratus 
and/or fog from cumulus type clouds. There was no apparent 


gray shade variation detectable between areas reporting fog 
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and those reporting code figure 6 for "G4" tise Stratus 

in a more or less continuous layer and/or ragged shreds, or 
both, but no stratus fractus of bad weather). feaeue 
Pierre wene discernible differences between areas of reported 
fog and other types of stratiform clouds where some vertical 
motion is implied, e.g., stratocumuliform clouds (code 
meure 4, 5, or 8 for "Ci"). Moreover, when using simultan- 
Somer oRVIS and SRIR-DAY imagery together, subtle variations 
gn both visual and IR frequently did occur in the vicinity 
of fog/no-fog boundaries. In good imagery, i.e., with good 
mieyeshede Variation and relatively undistorted through 
paotographic enlargement, the boundaries were more obvious. 
aoce tea) to the belief that the EPRF DRDC "Slice,“’ DMSP 
threshold processing, or any similar digital processing 
would be extremely useful in discerning marine fog. 

After looking at the twenty-one days of satellite 
imagery and documenting the results, skill was developed in 
correlating gray shade variations with verifying fog. As a 
result a set of seven qualitative guidelines were formulated 
for discerning areas where the satellite imagery indicated 
miter og Was Occurring. Phe giiidellines are offered as 
Appendix A. | 

The next step was to re-screen the satellite imagery and 
MiomCocuientatton compiled duying the first screening in 
PeactetOmseclect a sub-stUdy period. Considering as many 
factors as possible, the period selected was 241800 GMT 


through 270000 GMT July 1973. This period included a set of 
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three SRVIS/SRIR-DAY mosaics and two SRIR-NGT mosaics. The 
first six qualitative guidelines were then applied as objec- 
tively as possible to the eight enlarged mosaics and areas 
mamroo, no-fog,'' and “frontal clouds'' (i1.e., multilevel 
mromds) were established for the five primary synoptic study 
mmes, 1.€., 241800 GMT, 250600 GMT, 251800 GMT, 260600 GMT, 
and 261800 GMT. The secondary times (250000 GMT, 251200 GMT, 
260000 GMT, 261200 GMT, and 270000 GMT) were used wherever 
necessary for continuity. The areas diagnosed as "fog" were 
then verified against the synoptic reports and the statistics 
were compiled as shown in Table II. For the five study 
periods, the results showed that roughly 70% of the fog 

reports fell into the areas deemed, from the satellite imagery, 
to be fog areas. The skill of discerning fog as compared to 
chance (Panofsky and Brier, 1968) was positive for all cases 
Cmeeptetor the SRIR-NGI for the first night period, 250600 GMT, 
Gieenewas -0-06. ihe skill scores for each time period for 
Giemedses Of including and excluding the "frontal" data are 
also shown in Table II. 

Generally, the results show that one can discern areas 
where fog exists with some degree of positive skill based on 
the guidelines offered. In addition, it is believed that 
_Wimimictretmilmagery the score could have been even higher. 
Bemecnmenese nesults Ved to the use of the digital NOAA-2 
data, from which the mosaics were derived, for more objectively 


Giscernineg Marine fog and its horizontal extent. 





fee Ce ROGCESoING NOAA-Z DIGITAL COUNT DATA AND THE USE OF 
DIGITAL INFORMATION FOR DISCERNING MARINE FOG 


The NOAA-2 digital count data were provided on magnetic 
mame ald Coniigured as shown in Figure 2. The digital count 
meeaswere Obtained only for the sub-study period. The data 
fees, ranging from 0 to 255, were for a 2048 x 2048 grid- 
memieearray which coincided with the 20.7 x 20.7 cm 
digitally composited mosaics. A schematic of the mosaic 
background is provided in Figure 5 with the area of study 
within the shaded portion. The rectangular box incorporating 
mirem@Oceanic part of this area represents the outside bound- 
metecmot the data subpoints extracted from the 2048 x 2048 
array for detailed processing. The rectangular region is 
Oriented, for convenience of processing, so that the rows 
and columns of the sub-grid area are parallel to the rows 
and columns of the larger grid. The actual rows and columns 
Poieere sUp-orid area relative to the larger grid were 
eeoecmianed by Simple linear ratios involving 1) the linear 
eoeicesrrom the edges of the larger grid to the smaller 
Mekcmpomunidaries, 2) the total length of the larger grid 
(20.7 cm), and 3) the total number of grid spacings for the 
large grid (2047). For example, the uppermost grid row was 
Pome Dy weakingsthne ratio (9.3 cm):(20.7 cm)::y:(2047 spacings) 
MimmoolviIne Obey, in this case, y equals 920 grid spacings. 
(Vsti waseaetecrmined that the rows ranged from 921 to 1500 
and the columns from 446 to 792. This gave a 347 x 780 grid 


array for the sub-grid area. 


31 





The same geographical area covered by the sub-grid area 
was Outlined on a 1:15 million polar stereographic chart and 
Pienemedsured to determine its size and the expansion factor 
fiwy7so. the area, outlined on the 1:15M, measured 18 x 
10 3/4 inches. This rectangular area was then divided into 
a 1/4 inch meshlength grid. In order to transform the array 
of values from the finer mesh (approximately 1/32 inch at 
meeoMyeto the 1/4 inch meshlength grid at 1:15M, it was neces- 
‘wet ouextract roughly every eighth grid point. The 1/4 
inch meshlength gives a resolution of approximately 51 nmi 
at 60N and 45 nmi at S30N. Computer software was written to 
Pee all count values within the sub-area and also to print 
every eighth value, corresponding to the 1/4 inch meshlength 
pela, the 1/4 inch meshlength values were printed at 1/4 inch 
intervals so that clear acetate could be overlaid and a hand 
analysis made of the values. 

A rough hand analysis was then made of the 1/4 inch mesh 
@mieteal Count data for one SRVIS, SRIR-DAY, and SRIR-NGT 
iecalemto actermine a realistic range of values for computer 
analyzing the data at finer meshlengths. A range of 60 to 
MiGediettab counts con both SRVIS and SRIR was found to be 
mealistic. Attempts were then made to analyze the count data 
for the full 1/32 inch and 1/16 inch meshlength grids at 20- 
count intervals between 60 and 200. However, the analysis 
Soutware (NPS software-library routine "CONTUR") was not 
capabie of processing at these meshlengths. The smallest 


successfully analyzable meshlength was 1/8 inch (approximately 
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yom res@lucizon). It was tleen found that the 20-count interval 
my Ses Were basically too “noisy” for practical usage. 
icherere. thcmnexteseep was to correlate Metal count 
femmes with verifying fog reports. The range of count values_ 
for two SRVIS, two SRIR-DAY, and two SRIR-NGT cases (241800 GMT 
through 260600 GMT July 1973) are shown in Table ie re 
Maewalso tound that, for the SRVIS, digital counts less than 
Simeorresponded to generally cloud-free conditions; Cogan 
and Willand (1974) found that values less than 80 corresponded 
momvery tew or no clouds. 

Based on these values, the SRIR data were analyzed at 
wee seo oO. 120, 175, and 190 digital counts (roughly 224.0K, 
moromOnwweco' OK, and 298.0K, respectively). The SRVIS data 
Were analyzed at 60, 120, and 180 counts. Samples of the 
SRVIS, SRIR-DAY, and SRIR-NGT analyses are shown in Figures 
S$; 9, and 10. For the daytime data, ocean areas in which the 
SsRIR-DAY Count values were between 120 and 175 and which also 
Pomemiigerc (lcd -trec wae. elc ORVIS count values were greater 
than 60) were diagnosed as fog areas. For the nighttime data, 
pecrsmilene the okhPR-NGl count values merely fell between 120 
and 175 were diagnosed as fog areas. The verification scores 
Wem tnese § CUt-OLL digital Borne values are contained in 
aplenty ——ror salh bit ome time period, a considerable improve- 
jemiiesimitl Scores occurred; the first SRIR-NGI case even 
Seomed ad euancge Of Sign an the skill Scores (1.e., "chance" 
was beaten). Even more noticeable is the fact that all 


verifying fog reports fell into an area diagnosed as either 
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mateo trontal™ {(1.¢., multilevel clouds). More speci- 
mecdtivyy 3/2 Of the verifying fog reports fell into a fog- 
@eetonated area. Additionally, it was found that areas of 
multilevel clouds ("frontal") were more precisely defined in 
the computer produced display than in the enlarged mosaic 
imagery. 

iis. it) Was concluded that, for this particular sample 
of data, it was possible to diagnose marine fog with a good 
degree of positive skill using properly processed NOAA-2 


ioeital count data. 


ieee es -E ROCESSING DMSP SATELLITE IMAGERY 

tiene prelkiminary studies for this project, the DMSP 
satellite imagery in the form of transparencies, showed great 
MaecmlESses aS a tool for detecting lower tropospheric moisture 
omamstratus and fog. The transparencies, for which the data 
were earth-locatable (usually along the west coast of North 
America), were compared to verifying transient-ship synoptic 
weather reports. Generally the I-format IR imagery showed 
Smimebanre results as the NOAA-Z SRIR except for resolution. 
The H-format visual imagery and, particularly, the V-format 
Visual imagery showed some rather remarkable features which 
VoLecmiOusamparent an the NOAA-2 SRVJS. For example, along 
PicmperNpinery Of reported Stratus/fog bands, the white 
stratus/fog diminished to a very light shade of gray (lighter 
than the darker clear areas over the open ocean). This 


light shade of gray invariably extended around the entire 


34 





periphery of the stratus or fog and frequently extended down- 
wind of the cloud mass for several hundred kilometers. 

iitmeotiic cespcct mit 15 to be noted that the upper portion 
Peewee UMoP visual spectral interval in Figure 1 overlaps the 
water vapor region of the atmospheric absorption spectra. 
In the "positive" visual imagery, one would expect darkening 
of the gray shades if attenuation due to water vapor absorp- 
tion was dominating. However, the imagery actually showed 
a lightening of the gray shades where one would expect to have 
greater concentrations of lower tropospheric moisture. This 
meee to the belief that the attenuation (if it is in fact 
attenuation) is due to scattering, diffraction, or reflection 
from condensation nuclei, aerosols, or small cloud droplets. 
Meeordine to Rayleigh Theory, the ideal ratio of the diameter 
of a spherical particle to the wavelength of incident radia- 
Mmii@mrOcescattering 15 0.1. Thus for the DMSP visual spectral 
mocponse o1 0.4 to l.lum, particles of size of 0.04 to 0.liun 
would cause Rayleigh scattering. This would include such 
particles as Aitken nuclei, haze, and small condensation 
f(iceli@emeror two particular periods studied, there were veri- 
fying ship reports across the boundaries of the bright-white 
stratus/fog and the light-gray shaded area; in both cases 
ships beneath the brighter portion reported fog and ships 
beneath the light-gray areas reported haze. 

Thus, it is suspected that the DMSP satellite visual 
imagery may be extremely helpful in locating fog boundaries. 
Dpa@ikemoene NOMA-2 SRVIS digital data, it is quite possible 


Bisel cut-off processing value may be determined for the 


So 





DMSP digitally mapped 
the boundaries of fog 
then be used for more 
regions to be used in 


As of the writing 


Visual data as a means of isolating 

and no-fog. Hence, the visual data may 
than merely determining cloud-free 
Comjumetlon with IR data. 


of this thesis, the DMSP satellite 


data, which was processed by AFGWC, had not been received; 


therefore, a thorough 


inivestweatton into the DMSP satellite 


data for the entire area of study was not performed. 
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III. CONCLUSIONS AND RECOMMENDATIONS FOR FURTHER STUDIES 


It was the objective of this study to determine the 
feasibility of using meteorological satellite data for dis- 
Semmine Marine £08. It is believed that sufficient evidence 
is presented, from the small sample taken in this study, to 
justify further investigations into the real-time use of 
Satellite data for diagnosing marine fog. Both the qualita- 
tive guidelines in Appendix I and the numerical processing 
techniques used in this study showed the capability of 
mememesve Skill in diagnosing marine fog. Furthermore, it 1s 
believed that refinements to these techniques could 
assuredly improve the diagnostic skills. 

Therefore the following recommendations are offered for 


imeure studies: 


1) Develop computer software to earth-locate the 
digital count values for NOAA-2, given any latitude 
and longitude of verifying "ground truth" data; 
then, more accurately determine the appropriate 
ducttaie count values (SRVIS and SRIR) correspond- 
PrmLOuverityinc LO reports. Accumulate a large 
data sample of these values for the entire 21-day 
ivr perrtOdeanad Statistically derive the best 
€ut-o1f values for analyzing the digital count 


data. Determine the latitudinal variation, if 


oF 





any, of the count values corresponding to verify- 
Mig tog map iysethese ent-off values to the digi- 


PAmdelnon the £411 21-day period. 


2) Develop similar digital processing software 
for the DMSP mapped digital data whenever it is 
received from AFGWC. Also, test the qualitative 
guidelines in Appendix A and amend them as 


deemed necessary. 


3) Perform a real-time test of NOAA-series DRSR 
feta da live recording of the current NOAA 
Teco mcr Converting the data using the 
DRDC at FNWC, and digitally processing the data 
using the DRDC support software at EPRF (in parti- 


Heit zo stiee Stnce processing technique). 


Weeenerrormnisa real-time test of DMSP satellite 

data by requesting threshold processing of 

direct-readout DMSP data received at the Naval 

Lomticheockrvice racility, san Diego, California 

Homeric wcame Nertogds as the NOAA series data 

described in 3) above. 

If the numerical diagnosis of marine fog shows continued 

success, eventuaily a marine fog climatology of visual and 
IR digital count "cut-off" processing values can be compiled 


for every month of the year for various maritime regions of 





the world using the same or similar techniques performed or 
recommended in this study. Whenever this stage in develop- 
iemiels reached, a reasonable marine fog analysis may be a 


reality at last. 
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Peden ss oatiple OULpUt of Sorted, reconfigured, and categor- 
iMedwetransilent-Siip SyNoptiec-time reports contain- 
ing fog-related parameters. Output includes the 
iene londmemeadiorGalileGton (IREGS) of the ship, 
the fog category (CAT), the latitude (LAT) and 
honemeucemGONG mom the ship's position, the pres- 
ent weather code (WW), past weatrer code (W), visi- 
bility code (VV), total cloud amount/low-mid cloud 
amount (N/NH), and cloud types (LO, MID,HI) (X"s 
Mmiteaces tide the sky lseobscured ). 
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Code figure 
(ww=present weather 


Depth of 
fog at sea 
















Category Wee eerner) DeScrilpimaon 
a ww=10 with or without Light fog at station 
W=4 at time of observa- 
hon 
2 Ww=11, 12 with or Shallow fog at sta- 
without W=4 tion at time of 
observation 
3 ww=28, with or without 





W=4 






Fog at a distance 
(40) or at the sta- 
tion. (4i1=49) at 
time of observation 


ww=40-49 with or 
without W=4 






Fog or thick haze or 
smoke within a period 
1-6 hours prior to 
observation 






W=4 
ww? 10,11,12,28,40-49 






£ 
pe 
pa 
Fog at station not 
at time of observa- 33 
Elana bUbewa thin 
past hour 
S 


Table I. Categories used for cataloging transient-ships 
neoponcmic moo-relatcd Code-ficgures in synoptic- 
time weather reports. 
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135-176 98-190 ne 3: 
70-112 iso =220 ie 0 elie 















Fog beneath 
frontal bands 
or multilevel 
clouds 






Table III. Range of NOAA-2 SRIR and SRVIS 
digital counts associated with 
Verio marines fogs tor the 
period 241800 GMT through 
260600 GMT July 1973. 
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APPENDIX A 


QUALITATIVE GUIDELINES FOR OUTLINING AREAS OF 
MARINE FOG USING NOAA SRVIS, SRIR-DAY, AND 
SRIR-NGT IMAGERY IN CONJUNCTION WITH 
Bree NoLENT-SHiP SYNOPTIC DATA AND SEA-LEVEL PRESSURE ANALYSES 


Outline areas in the SRVIS and SRIR-NGT which appear 
eloud-free (should appear dark). 

iiiime areas in SRVIS where the cloud texture appears 
“"bumpy'' (which implies cumuliform or stratocumuliform 
clouds) and note the areas where clouds have a smooth 
Peace im SkIR-NGI, the bumpy-textured areas should 
appear as patches of alternating light and dark gray 
cides.  sthie smooth textured clouds will appear as a 
nearly continuous light gray shade. 

Sitlime, as well as possible, the total cloud band of 
frontal zones in SRVIS or SRIR-NGT (including areas 
with apparently multilevel clouds). 

Observe the SRIR-DAY imagery in the areas declared as 
Weloud=e:tree’’ in Step I. Note the gray shade of this 
reference "'clear-column" (for SRIR-NGT, this should 
have been done in Step 1). 

Note the areas in the SRIR imagery which display smooth 
textured clouds (excluding frontal bands) that radiate 
nimmecdark fo medium cray shade (i15e., slightly lighter 
Gimmes clear column), and outline this area as an 


aEeamOormtog. ff a white appearance occurs in the 
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medium gray shaded area, assign “no-fog" to the region 
under that white area. 

Observe the portions of the frontal band along the 
leading and trailing edges. If there are regions with 
emeceichtly darker pray shade than the cold, multi- 
Mayenca Clouds of the main frontal band and slightly 
lighter than the clear-air masses in advance of and 
Memenad the tromtal band, outline these areas as areas 
Gi LOe. 

SBenpare the areas designated as “fog" and "no-fog" 
with available transient-ship synoptic data and the 
euaneme sea-level pressure analysis. Adjust the 
satellite-designated fog areas to be compatible with 


circulation/thermal features and marine observations. 
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